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Reliable preparation of entanglement between 
distant systems is an outstanding problem in 
quantum information science and quantum com- 
munication. In practice, this has to be accom- 
plished via noisy channels (such as optical fibers) 
that generally result in exponential attenuation 
of quantum signals at large distances. A spe- 
cial class of quantum error correction protocols — 
quantum repeater protocols — can be used to 
overcome such losses. In this work, we intro- 
duce a method for systematically optimizing ex- 
isting protocols and developing new, more effi- 
cient protocols. Our approach makes use of a dy- 
namic programming-based searching algorithm, 
the complexity of which scales only polynomially 
with the communication distance, letting us effi- 
ciently determine near-optimal solutions. We find 
significant improvements in both the speed and 
the final state fidelity for preparing long distance 
entangled states. 



INTRODUCTION 

Sequential decision making in probabilistic systems is a 
widely studied subject in the field of economics, manage- 
ment science and engineering. Applications range from 
problems in scheduling and asset management, to control 
and estimation of dynamical systems [ ]. In this paper 
we make the first use of these techniques for solving a 
class of decision making problems that arise in quantum 
information science [2, •'->]. Specifically we consider the op- 
timal design of a so-called quantum repeater for quantum 
communication. Such repeaters have potential applica- 
tion in quantum communication protocols for cryptog- 
raphy [4, 5, 6] and information processing [7], where en- 
tangled quantum systems located at distant locations are 
used as a fundamental resource. In principle this entan- 
glement can be generated by sending a pair of entangled 
photons through optical fibers. However, in the presence 
of attenuation, the success probability for preparing a 
distant entangled pair decreases exponentially with dis- 
tance [8], 

Quantum repeaters can reduce such exponential scal- 
ing to polynomial scaling with distance, and thus provide 
an avenue to long distance quantum communication even 
with fiber attenuation. The underlying idea of quantum 



repeater [9, 10] is to generate a backbone of entangled 
pairs over much shorter distances, store them in a set 
of distributed nodes, and perform a sequence of quan- 
tum operations which only succeed with finite probabil- 
ity. Purification operations [11, 12] improve the fidelity 
of the entanglement in the backbone, while connection 
operations join two shorter distance entangled pairs of 
the backbone to form a single, longer distance entangled 
pair. By relying on a quantum memory at each node 
to let different sections of the repeater re-attempt failed 
operations independently, a high fidelity entangled state 
between two remote quantum systems can be produced 
in polynomial time. A quantum repeater protocol is a 
set of rules that determine the choice and ordering of 
operations based upon previous results. An optimal pro- 
tocol is one that produces entangled pairs of a desired 
fidelity in minimum time within the physical constraints 
of a chosen implementation. 

The complexity of finding the optimal repeater pro- 
tocols can be understood via the following analogous 
example problem [1]: given a sequence of rectangular 
matrices M\Mi . . . M n , such that is dk x rffc+i di- 
mensional, find the optimal order of multiplying the ma- 
trices such that the number of scalar multiplications is 
minimized. This is a typical example of a nesting prob- 
lem, in which the order in which operations are carried 
out effects the efficiency. For example, if Mi = 1 x 10, 
M 2 = 10 x 1 and M 3 = 1 x 10, then (M 1 M 2 )M 3 takes 
only 20 scalar operations, while Mi(M2M 3 ), requires 
200 scalar multiplications. A brute force enumeration 
of all possible nesting strategics and evaluation of their 
performance is exponential in n. To solve this prob- 
lem more efficiently, we observe that the optimal nesting 
strategy (Mi ...(...)... M p )(M p+ i . . . M„) should carry 
out the solution to its subparts optimally, i.e., the nest- 
ing (Mi ...(...)... M p ) should represent the best nest- 
ing strategy for multiplying M1M2 . . . M p . This is the 
well-known dynamic programming strategy [ ], in which 
one seeks to optimize a problem by comparing different, 
already optimized sub-parts of the problem. Dynamic 
programming enables us to find the optimal solution to 
the original problem in time that is polynomial in n. 

Quantum repeaters also have a nested (self-similar) 
structure, in which shorter distance entanglement is used 
to create longer distance entanglement, which is then 
used in turn for further extending the distance between 
entangled pairs. This structure allows us to use the 
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FIG. 1: Quantum repeater scheme from Refs. [9, 10] 
(BDCZ scheme), (a) In a typical realization with N + 
1 = 9 nodes, the number of qubits per nodes is bounded by 
21og 2 2JV = 8. (b)-(d) Two entangled pairs with distance 
1 are connected (orange rounded rectangle) at node #1 to 
produce an entangled state with distance 2, which is stored 
(purple arrows) in the qubits at higher level, (e)-(g) Another 
entangled state with distance 2 is produced to purify (purple 
arrows) the entangled state stored in qubits at higher level. 
Similarly, entangled states with distance 2™ can be connected 
to produce entangled state with distance 2 n+1 , which may be 
further purified, as indicated in (a). 

methods of dynamic programming to find optimal nest- 
ing strategies for designing quantum repeater protocols. 

We now proceed to detail the specific optimization 
problem, then discuss our dynamic programming solu- 
tion to the problem. We next examine two representa- 
tive schemes that we wish to optimize [the Briegel et al. 
scheme (BDCZ scheme), Refs. [9, 10], and the Childress 
et al. scheme (CTSL scheme), Refs. [13, 14]], and find 
significant improvements in both preparation time and 
final fidelity of long distance entangled pairs. 

DYNAMIC PROGRAMMING APPROACH 

General Quantum Repeater Protocol. Quantum 
repeater protocols have a self-similar structure, where the 
underlying operations at each stage of the repeater have 
the same basic algorithms. In other words, the structure 
of the problem remains the same at each stage, while 
the parameters can be different. A generic quantum re- 
peater consists of three kinds of operations: entangle- 
ment generation, entanglement connection, and entan- 
glement purification. Entangled pairs are first generated 
and stored over a short distance Lq. At the first nesting 
level, two entangled pairs of distance Lq can be extended 



(a) 



o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 





1 


2 


k 


k+1 


n-2 


n-1 


n 




(c) 

O — — K> O j O 

1 k~~k+l n-1 n 

(d) ^ 

o o o "b 

(e) ^ 7 -_- 

o" O — — K> *j O O 

1 l" "l+l n-1 n 

, 

' v O-*- o^ 'o -►"oj 

(g) 

o e — o o 

(f) 

o (0^70+-. -.— ^oh^o/ ~b 

12 n-2 n-1 n 
(f") 



FIG. 2: Quantum repeater scheme from Refs. [13, 14] 
(CTSL scheme), (a) This scheme has exactly two qubits 
per node. The communication qubits (green nodes) are used 
for entanglement generation and short-term storage; the stor- 
age qubits (black nodes) are used for long-term storage, (b) 
With the help of local gates (orange solid rounded rectan- 
gles) between communication and storage qubits, the entan- 
gled state between communication qubits can be used to im- 
plement entangling gates (e.g. the Controlled NOT gate) 
between storage qubits from neighboring nodes. The effec- 
tive remote gate is highlighted by the orange dashed rounded 
rectangle. Such remote gate is sufficient for entanglement 
connection and purification of storage qubits. The commu- 
nication qubits are omitted in the following plots, (c)-(d) 
Entanglement connection to produce an unpurified entangled 
state with distance n. (e)-(g) Entanglement connection to 
produce an unpurified entangled state with distance n — 2 to 
purify the entangled state with distance n. (f')(f") Illustra- 
tion of multi-level pumping. An entangled state with distance 
n — 4 is used to purify an entangled state with distance n — 2, 
and the (purified) latter is used to purify an entangled state 
with distance n. The only difference between (f) and (f) is 
the fidelity of the entangled state with distance n — 2. The 
latter has higher fidelity. 
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to distance L\ ~ 2Lq via entanglement connection [6]. 
Due to limited fidelity of the short pairs and the im- 
perfections from the connection operations, the fidelity 
of the longer pair produced by connection is generally 
lower than those of the shorter ones. Nevertheless, the 
fidelity of the longer pair can be improved via entangle- 
ment purification, which is able to extract high-fidelity 
entangled pairs from an ensemble of low-fidelity ones us- 
ing operations that are local (restricted to qubits within 
a given node) [11, 12]. An efficient approach of entangle- 
ment purification is entanglement pumping [9, 10], which 
purifies one and the same pair by using low-fidelity pairs 
with constant fidelity. [ ] Thus, at the (fc + l)th nest- 
ing level, the three underlying operations (preparation at 
distance L^, connection, and purification) lead to prepa- 
ration at a distance Lk+i ~ 2Lk- [30] 

Inductive Optimization. We now define the opti- 
mization problem: 

Def: For given physical resources, desired distance 
ifinah an d final fidelity Ffi na i, an optimal protocol 
minimizes the expected time to have an entangled 
pair of fidelity F > Fft na \ at a distance L > Lfi na i- 

To solve this optimization problem, the choice of pa- 
rameters for the quantum operation cannot be viewed in 
isolation — one has to trade off the desire of low present 
cost (in terms of time) with the undesirability of high fu- 
ture costs. If one tries to enumerate and test all possible 
adjustable parameters, the complexity to search for the 
optimized implementation scales at least exponentially 
with the number of repeater nodes. A simple example 
is provided if we make our only adjustable parameter 
choosing between zero and one purification step at each 
stage of the protocol. For the BDCZ scheme with 128+1 
repeater nodes, there are already 2 128 > 10 38 possibili- 
ties, which is beyond the capability of current computers. 
Thus a systematic searching method is needed to find the 
optimized implementation out of such a huge parameter 
space. 

Based on the above self-similar structure, we may 
express the optimized protocol to produce long entan- 
gled pairs in terms of a set of optimized protocols for 
producing shorter pairs. The general searching proce- 
dure can be performed inductively, as detailed in Table 
I. We make a discrete set of target fidelities (see Sup- 
porting Methods for details), F = {/i, fo, . , . , f q }, such 
that only a finite number of different optimal protocols 
with shorter distances need to be developed. The com- 
plexity for each step of our optimization procedure is 
shown in the table; the full procedure scales as 0(q 2 n 2 ), 
where the 0(n) repetitions of step 3 take the most 
time. [ ] In practice, we found the full search of step 
3 to be unnecessary — the search can be restricted to 
pairs of distance n/2 ± C(log(n)), leading to complex- 
ity 0(q 2 nlog(n)). 



1. Find and store implementations that optimize the av- 
erage time (for all fidelities fi, ■ ■ ■ ,fq) with dist = n = 1, 
taking O(q) 

2. Assume known optimized implementations (for all fideli- 
ties) with dist < n 

3. Find optimized implementations to produce unpurified 
pairs (for all fidelities) with dist — n + 1 by trying (connect- 
ing) all combinations of known optimized implementations 
with dist < n, with complexity of order 0(q 2 n) 

4. Find optimized implementations to produce purified 
pairs (for all fidelities) with dist = n + 1 by trying all com- 
binations of unpurified pairs with dist = n + 1, pumping for 
m = 0, 1, 2, 3, • • • times; complexity goes as C(m max g 2 ) 

5. Store the optimized implementations (for all fidelities) 
with dist = n + 1, based on step 4. 

6. Replace n by n + 1, and go to step 2. 



TABLE I: Inductive search using dynamic programming 

Repeater Schemes and Physical Parameters. 

So far, we have only taken a general perspective in ex- 
plaining quantum repeater protocols and describing the 
procedure of inductive searching using dynamic program- 
ming. In this subsection, we specify the parameters to 
be optimized by examining the schemes of the quantum 
repeater, physical restrictions on entanglement genera- 
tion for current techniques, and the error models of lo- 
cal quantum gates. Only with a functional relationship 
between physically adjustable parameters and repeater 
operation outputs, can we find the optimized implemen- 
tations for procedure 1, 3 and 4 in Table I. 

There are several different schemes for building a quan- 
tum repeater that differ primarily in the amount of physi- 
cal resources utilized. For example, in the BDCZ scheme 
[9, 10] (Fig. 1), the maximum number of qubits in the 
quantum memory (to store intermediate states for con- 
nection and purification) required for each repeater node 
increases logarithmically with the total number of re- 
peater nodes. In the CTSL scheme [13, 14] (Fig. 2), 
an efficient way to use quantum memory is proposed, 
and only two qubits are needed for each node, regard- 
less of the total number of repeater nodes. One of the 
two qubits is called the communication qubit, which is 
optically active and can be used to generate entangle- 
ment with other communication qubits from neighboring 
nodes. The other qubit is called the storage qubit, which 
can be used to store quantum state over very long time. 
As shown in Fig. 2(b), with the help of local gates (or- 
ange solid rounded rectangles) between communication 
and storage qubits, the entangled state between commu- 
nication qubits can be used to implement tclcportation- 
based gates (e.g., the Controlled NOT gate) between 
storage qubits from neighboring nodes [7, 24]. Such re- 
mote gates (orange dashed rounded rectangle) are suffi- 
cient for entanglement connection and purification of the 
storage qubits; communication qubits are providing the 
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necessary resource mediating the gates between remote 
storage qubits. For clarity, we will omit the communi- 
cation qubits in the following discussion, but still keep 
track of the mediated operation between remote storage 
qubits. 

To model errors in the physical operations, we need 
to introduce a number of parameters determined by the 
quantum hardware. For entanglement generation, the 
relationship between the fidelity of the elementary pair, 
Fq, and the generation time, r e , depends on the physi- 
cal parameters (such as the signal propagation speed, c, 
the fiber attenuation length, L aU , the efficiency of single 
photon collection and detection, e, and the distance of el- 
ementary pair, Lq) and the specific approach to generate 
entanglement. For example, for the entanglement gener- 
ation approach using scattering as proposed in Refs. [13, 

, , , f r r r ,r l 2 ^)/^ 

], F = F (r e ) = 



1 



1 



For entanglement connection and pumping, the domi- 
nant imperfections are errors from measurement and lo- 
cal two-qubit gate, which we model with a depolarizing 
channel. In particular, the model for measurement is 
quantified by a reliability parameter [9, 10], rj, which is 
the probability of faithful measurement. For example, a 
projective measurement of state |0) would be 

P = r?|0) (0| + (l-r?)|l)<l|. 

Similarly, the model for local two-qubit gate is charac- 
terized by a reliability parameter [9, 10], p. With prob- 
ability p, the correct operation is performed; otherwise 
the state of these two qubits is replaced by the identity 
matrix [9, 10]. For example, the action on a two qubit 
operation [Zy would be 



l-p 



Try [p] 



where Tr^ [p] is the partial trace over the subsystem i 
and j, and Iij is the identity operator for subsystem i 
and j. Generally, the reliability parameters (77 and p) 
should be reasonably high (i.e., above some thresholds 
[9, 10]), so that the suppression of error from entan- 
glement pumping dominates the new errors introduced 
by entanglement connection and entanglement pumping. 
[32] 

Optimization Parameters. We now list the ad- 
justable parameters we can optimize over during proce- 
dures 1, 3 and 4 in Table I. 

1. During the entanglement generation, there is free- 
dom to choose the generation time r e , which is de- 
termined by the success probability and the com- 
munication time. Generally, the higher the success 
probability, the shorter the generation time and the 
lower the fidelity of the entangled state, so the gen- 
eration time and the fidelity should be balanced 
[13, 14]. 



2. During the entanglement connection, the distances 
of two shorter pairs can be adjusted, while the total 
distance is kept unchanged. 

3. During entanglement purification, the number of 
steps is also adjustable, which should balance the 
gain in fidelity and the overhead in time. 

Additional Operations. Besides the above opera- 
tions from the original quantum repeater schemes, there 
are some additional operations that might be useful. 
For example, we may skip several intermediate repeater 
nodes (node- skipping) to generate entanglement between 
distant nodes directly with a substantially lower suc- 
cess probability. Also, during entanglement pumping, we 
might consider multi-level pumping [15], which is to nest 
several levels of entanglement pumping together before 
the next level of entanglement connection (Fig. 2(f)). 
Multi-level pumping can produce entangled pair with 
higher fidelity than single-level pumping. Such additional 
operations can be easily incorporated into the search pro- 
cedures 1, 3, and 4 in Table I. We will show that the dy- 
namic programming approach can use these additional 
operations appropriately, to reduce the average time, ex- 
tend the upper bounds for achievable final fidelity, and 
even improve the threshold for the reliability parameters 
of p and 77. 



RESULTS AND DISCUSSION 

Improvement of BDCZ and CTSL Schemes. 

With procedure as listed in Table I, we implemented a 
computer program to examine the mean time to prepare 
entangled pairs and to search according to our dynamic 
programming prescription through the parameter space 
outlined above. We looked for optimal protocols for a 
quantum repeater for all distances < 1280 km and tar- 
get fidelities > 0.8. Unless otherwise specified, we use 
L at t — 20km, e — 0.2, and 77 = p = 0.995 for the rest of 
the discussion. We first fix Lq = 10km, and we will 
consider the optimization of Lq to justify such choice 
later. To visualize the results, the profile of the opti- 
mized time (smooth surface) is plotted in Fig. 3(a) (b) 
with respect to the final distance (from 10km to 1280fcm) 
and the fidelity (from 0.90 to 0.99) for both the BDCZ 
and the CTSL schemes. The calculation optimizes over 
the elementary pair generation (both distance and gen- 
eration time), the connecting positions, and the number 
of pumping steps, with spacing between neighboring re- 
peater nodes of 10km; both additional operations (node- 
skipping and multi-level pumping) are also included for 
the optimization. For comparison, the unoptimized time 
profiles (meshes) for the BDCZ and the CTSL schemes 
are also plotted. The unoptimized protocol assumes fixed 
elementary pair fidelity (Fq = 0.96 and 0.99 for BDCZ 
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FIG. 3: Plots of time profiles and improvement fac- 
tors. Speed-up in time associated with various final distance 
and fidelity, (a) t (Ffi na i, L) for unoptimized (meshes) and 
optimized (smooth surface) implementations of the BDCZ 
scheme; (b) for the CTSL scheme, (c) Pseudocolor plot of 
the improvement factor, t unop t/t op t, for the BDCZ scheme; 
(d) for the CTSL scheme, in the region (Ff ina i > 97.5), the 
improvement factor tunopt/topt — > oo. The default parameters 
are L a tt ~ 20km, e = 0.2, and p — r\ — 0.995. 

and CTSL, respectively), simple connection patterns (de- 
tailed in Ref. [9, 10] and [13, 14]), and constant number 
of pumping steps. 

As expected, the unoptimized protocol always takes 
longer time than the optimized protocol for the same fi- 
nal distance and target fidelity. Time profiles for the un- 
optimized protocols have stair-like jumps in Fig. 3(a) (b). 
For the BDCZ scheme (Fig. 3(a)). The jumps occurring 
with increasing distance (occurring at distances L/Lq = 
2 P + 1 = 1,3,5,9, 17, 33, • • • ) are the results of time over- 
head from the additional level of connection; the jump 
occurring with at Ffi na i = 0.98 is due to the sudden 
change in the number of pumping steps from 1 to 2. Sim- 
ilarly, for the CTSL scheme (Fig. 3(b)), the two jumps 
are due to the change of the number of pumping steps 
from 1 to 2 and finally to 3. For the optimized protocols, 
the time increases smoothly with increasing final distance 
and fidelity 

The improvement factor (i.e., the ratio between the 
times for unoptimized and optimized protocols) is plot- 
ted for both the BDCZ and the CTSL schemes in 
Fig. 3(c) (d). As we might expect, the previously men- 
tioned jumps lead to sharp stripes where the improve- 
ment factor changes discontinuously. There are several 
regions where the optimization gives significant improve- 
ment. For example, for the BDCZ scheme, the vertical 
bright stripes indicate that the optimization provides a 
time-efficient way to generate entangled pairs for distance 
(2 P + <5 + ) Lq (with S + > 0), gaining a factor of about 10; 



FIG. 4: Plots of time profiles and improvement fac- 
tors. The subplots are arranged in the same way as Fig. 3. 
Local operations have lower reliability parameters, p = n = 
0.990. (a)(c) For the BDCZ scheme, the optimization proce- 
dure only improves the speed of the quantum repeater, and 
does not extend the achievable region in the F-L plot. (b)(d) 
For the CTSL scheme, for distances longer than 200fcm, the 
improvement factor, t unopt /t op t — > oo. Here, the reliability 
parameter (p = r\ — 0.990) is insufficient to create distant en- 
tangled pairs with the unoptimized implementation, but the 
optimized implementation (with multi-level pumping) is still 
able to create high-fidelity distant entangled pairs, because 
multi-level pumping lowers the threshold of the reliability pa- 
rameters for the CTSL scheme. 

the horizontal bright stripes indicate that efficiently ar- 
ranging the number of pumping steps can also speed up 
the scheme by a factor of about 30 or even more. For 
most of the optimized protocols, a distant pair is divided 
into two shorter pairs with similar distance and fidelity 
(symmetric partition), but occasional asymmetric parti- 
tioning can further reduce the time by about 10%. 

For the BDCZ scheme, the correspondence between 
jumps and stripes essentially accounts for all the fea- 
tures of the improvement plot (Fig. 3(c)). For CTSL 
scheme, however, besides the stripes, there is also a re- 
gion (with distance L > 100km and fidelity F ^ 97.5%) 
where the improvement factor is infinity — optimization 
not only boosts the speed, but also extends the upper 
bound of achievable fidelity for distant pairs. 

We also study the improvement for other choices of 
reliability parameters, p and 77, especially those values 
close to the threshold [9, 10]. Suppose the reliability pa- 
rameters are p = rj = 0.990. In Fig. 4(a) (c), we plot 
the speed-up in time associated with various final dis- 
tance and fidelity for the BDCZ scheme. For both (op- 
timized and unoptimized) protocols, the highest achiev- 
able fidelity is approximately 97.5% (compared to 99% in 
Fig. 3(c)), limited by errors from local operations. The 
improvement factor ranges between [1.5,600] (compared 
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FIG. 5: Examples. Two implementations with targeting final distance L = HLo and fidelity Ffi w i = 0.976, using the 
CTSL scheme. Each storage qubit is represented by "*". All the relevant entangled states are shown. The order to produce 
these entangled states are from bottom to top; states on the same row can be produced simultaneously. There are two kinds of 
entangled states - purified entangled states (type- A, solid black line) and unpurified entangled states (type-B, dashed blue line). 
On the left side of each purified entangled state, there is a label "Afc", and this number k indicates that this purified entangled 
states is obtained from k steps of entanglement pumping, (a) The unoptimized (left) implementation has three pumping steps 
after each entanglement connection, with average time of about Usee to produce the pair wanted, (b) The optimized (right) 
implementation is from optimization over pair generation time, connection position, and the number of pumping steps. The 
optimized choice of connection position does not necessarily break the long pair into two almost identical shorter pairs; for 
example, the entangled state pointed by the black arrow in the 9th row is obtained by connecting two very different shorter 
pairs in the row below. In addition, the possibility of multi-level pumping is also taken into account during the dynamic 
programming. As pointed by the red arrows, the pair of storage qubits in the third row pumps the pair in the second row, and 
the latter pumps the pair in the first row. The average time is about 1.5 sec for the optimized implementation, about 8 times 
faster than the unoptimized one. 



to [1,100] in Fig. 3(c)). Apart from these differences, 
the key features (horizontal and vertical stripes) of im- 
provement from optimization are very similar between 
Fig. 3(c) and Fig. 4(c). 

For the CTSL scheme, however, unoptimized and op- 
timized protocols behave very differently, when p = r\ = 
0.990. As shown in Fig. 4(b) (d), the unoptimized proto- 
col cannot effectively create entangled pairs for distances 
longer than 200/cm, while the optimized protocol is still 
able to efficiently create distant entangled pairs with very 
high fidelity. Thus our optimization lowers the threshold 
requirement for the CTSL scheme of quantum repeater. 

To understand the reason for the improvement of the 
highest achievable fidelity (Fig. 3(d)) and the parameter 
threshold (Fig. 4(d)), we examine the optimized protocol 
of CTSL scheme in the next two subsections. 

Comparison between Optimized and Unopti- 
mized Protocols. We first compare the detailed pro- 
cedures between two (optimized and unoptimized) pro- 
tocols using CTSL scheme to produce a pair with final 
distance L = HLq and fidelity Ffi na i = 97.6%, with de- 
fault reliability parameters p = rj = 0.995. We choose 



the highest fidelity achievable by the unoptimized pro- 
tocol, so that we will see almost all features that give 
improvements. The results for the unoptimized protocol 
(Fig. 5(a)) follows Refs. [13, 14] exactly, while the opti- 
mized one (Fig. 5(b)) is from our systematic search using 
dynamic programming. They differ in the following as- 
pects: (1) during entanglement generation, the optimized 
implementation generates elementary pairs with fidelity 
lower than 0.99 to reduce the generation time, and uses 
entanglement pumping afterwards to compensate the fi- 
delity loss; (2) during entanglement connection, the rule 
of producing long pair from two almost identical shorter 
pairs is slightly modified (e.g., the pair pointed by the 
black arrow in the 9th row is from connection of two 
quite different pairs in the 10th row); (3) the number of 
pumping steps after each connection varies from to 3 
for optimized implementation; (4) finally, the optimized 
implementation uses multi-level pumping, which will be 
discussed in detail in the next subsection. For clarity, the 
additional operation of node-skipping is suppressed in the 
optimization here. The overall average time is reduced 
from Usee to 1.5 sec, improved by a factor of 8. 
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FIG. 6: Results of Monte Carlo simulation. Monte 
Carlo simulation for unoptimized/optimized implementations 
for (a) the BDCZ scheme and (b) the CTSL scheme, with fi- 
nal distance 1280fcm and fidelity 0.97. The time distributions 
for distant pairs are plotted, with red (blue) bars for opti- 
mized (unoptimized) implementation. In each plot, the red 
(blue) diamond indicates the estimated time from average- 
time approximation for optimized (unoptimized) implemen- 
tation. The average-time approximation provides a good es- 
timate up to some overall factor (2 ~ 3) which is not very 
sensitive to the implementation. 



Note that our optimization results based on average 
time approximation (see Fig. 6 and Supporting Meth- 
ods) are confirmed by the Monte-Carlo simulation of the 
optimized protocols, verifying the substantial speed-up. 

Multi-Level Pumping. We now consider the ad- 
ditional operation of multi-level pumping in more de- 
tail. We discuss multi-level pumping only for the CTSL 
scheme, but not for the BDCZ scheme. (In the BDCZ 
scheme, to introduce multi-level pumping requires addi- 
tional quantum memory qubits.) In the original unopti- 
mized protocol [13, 14], the purified entangled state with 
distance n (between the Oth and the nth nodes (n > 5)) 
is produced by entanglement pumping, and the entangled 
states used for pumping (called pumping-pairs) are un- 
purified entangled states with distance n — 2 (Fig. 2(f)). 
The fidelity of these pumping-pairs with distance n— 2 are 
limited by the connection operation, which imposes an 
upper-bound for the fidelity of the purified pair with dis- 
tance n. The underlying restriction is that the pumping- 
pair is unpurified. 

We may lift this restriction by allowing the use of a 
purified pumping pair. This is multi-level pumping. For 
example, the pumping-pair with distance n — 2 may also 
be produced by entanglement pumping from pumping- 
pairs with distance n — 4 (Fig. 2(f)), and so on. By 
doing multi-level pumping, the fidelity of the pumping- 
pair with distance n — 2 is increased (Fig. 2(f)), and the 



same for the fidelity upper-bound for the entangled state 
with distance n. While multi-level pumping can increase 
the fidelity, it also slows down the repeater scheme. 

When the reliability of local operations is above the 
threshold for the unoptimized protocol (e.g. p = rj = 
0.995), we find that multi- level pumping is necessary only 
for the last two or three levels to the high-fidelity pair 
we want to produce. Such multi-level pumping can be 
identified in the optimized implementation - for example, 
as indicated by red arrows in Fig. 5(b), the pair of storage 
qubits in the third row pumps the pair in the second 
row, and the latter pumps the pair in the first row. On 
the other hand, when the reliability of local operations 
is below such threshold (e.g. p = rj = 0.990), multi- 
level pumping is needed almost after every entanglement 
connection. 

If we exclude the possibility of multi-level pumping 
in dynamic programming, the infinite improvement fac- 
tor for pairs with distance L > lOOfcm and fidelity 
F 97.5% in Fig. 3(d) would disappear. Similarly, in 
Fig. 4(d), without multi-level pumping, there would be 
no improvement of the parameter threshold, and even 
the optimized protocol could not efficiently create distant 
(L > 200km) entangled pairs. For the CTSL scheme, 
multi-level pumping not only enables us to prepare en- 
tangled pairs with very high fidelity, but also lowers the 
required threshold of the reliability parameters {p and rj) 
for local operations. Therefore, the flexibility to include 
additional operations in our dynamic programming pro- 
vides a new perspective on the optimization of quantum 
repeater schemes. 

Other Improvements. In addition to the previ- 
ously discussed features in the plots of improvement fac- 
tor, there is an overall improvement for all final distances 
and fidelities. Such overall improvement comes from the 
optimized choice of the distance (by node-skipping) and 
the generation time for each elementary pair used. Such 
overall improvement is about 1.5 (or 2^3) for the 
BDCZ (or CTSL) scheme, which indicates that the orig- 
inal choice of uniform distance Lq = lOfcm and initial 
fidelities Fq = 96% (or 99%) are quite close to the opti- 
mal. 

Finally, we consider if it is possible to gain some addi- 
tional speed-up if we are allowed to choose the location 
of the nodes of the quantum repeater. In order to an- 
swer this question, we discretize the distance into smaller 
units, e.g. lkm <C L att . Since the distance of each el- 
ementary pair is determined by the dynamic program- 
ming, the optimized location of the nodes can be inferred 
from the distances of the elementary pairs. We find that 
the speed-up due to optimization over the location of 
the nodes is fairly small, no more than 15% in time (for 
cases with final distances larger than 200fcm). Generally 
we find that as long as the node spacing is less than the 
attenuation length (Lq < L at t), a quantum repeater can 
be implemented almost optimally. 
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Experimental Implications. Throughout our 
analysis we have assumed relatively high fidelity of 
local measurements and operations (n = p = 0.995 or 
0.99) and memory times exceeding total communication 
times. Recent experiments with tapped ions [16, 17], 
neutral atoms [18], and solid state qubits [19] are already 
approaching these values of fidelity and memory times. 
At the same time, high initial entanglement fidelity 
(F w 96% or 99%) is also needed for the optimized 
protocols. Entanglement fidelity of about 90% can 
be inferred from recent experiments with two ions in 
independent traps [20]. While optimization procedure 
can yield protocols compatible with fairly low initial 
fidelity and high local error rates, in practice these errors 
introduce a large overhead in communication time. 

Besides the schemes considered here, there exist other 
quantum repeater schemes, in particular the Duan et al. 
scheme (DLCZ scheme) [21] that requires a smaller set 
of quantum operation and relatively modest physical re- 
sources. The original DLCZ scheme does not use active 
entanglement purification and hence cannot correct arbi- 
trary errors. In such a case, optimization is straightfor- 
ward and has been discussed in [21]. Recently, the DLCZ 
scheme has been extended to include active entanglement 
purification in order to suppress e.g. phase noises [22, 23]. 
The extended DLCZ scheme becomes very similar to the 
BDCZ scheme in terms of the self-similar structure. The 
technique of dynamic programming can be applied to op- 
timize the extended DLCZ scheme as well. 



CONCLUSION AND OUTLOOK 

We have demonstrated how dynamic programming can 
be a powerful tool for systematically studying the opti- 
mization of quantum repeater protocols. We find sub- 
stantial improvements to two specific repeater schemes [ , 
10, 13, 14]. Beyond searching for optimal choices in previ- 
ously known elements of the schemes (entanglement gen- 
eration, connection, and pumping), our systematic study 
can also incorporate more sophisticated additional oper- 
ations, such as node-skipping, multi-level pumping, and 
the flexible location of repeater stations. In particular, 
our multi-level pumping procedure extends the maximum 
achievable fidelity for distant pairs. It should be possi- 
ble to include additional possibilities to the optimization 
problem of quantum repeater, such as different choices of 
entanglement generation and possibly more efficient us- 
age of local qubits [25, 26]. It would also be interesting 
to study the optimization problem of quantum repeater 
with finite storage time of the quantum memory [27, 28]. 
Even the optimized protocols have a rather limited speed 
(corresponding to generation of one high-fidelity pair over 
1280 km in 1 ~ 100 s (see Fig. 6). Therefore, improve- 
ment of experimental techniques (to obtain higher lo- 
cal operation fidelities and more efficient atom-photon 



coupling) as well as development of new theoretical ap- 
proaches to speed-up quantum repeaters still remain an 
outstanding goal. Furthermore, the dynamic program- 
ming techniques may find application in other outstand- 
ing problems in quantum information science, such as the 
optimization of quantum error correction for fault toler- 
ant quantum computation. In particular, the optimiza- 
tion of the network-based quantum computation scheme 
with minimal resources [24] might be possible. 
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SUPPORTING METHODS: SHAPE 
PARAMETER APPROXIMATION AND 
AVERAGE TIME APPROXIMATION 

We use two important approximations throughout the 
analysis: the shape parameter approximation and average 
time approximation. 

In the shape parameter approximation [10, 13, 14], we 
use two numbers, the fidelity and the shape parameter 
(F, v) to classify non- ideal entangled states during our 
inductive optimization procedure. For a pair of entan- 
gled qubits, we assume that the density matrix is diag- 
onalized in the Bell basis p — diag (F\, F2, F3, F4), with 
F\ > F 2 > F3 > F4 by an appropriate ordering of the 
Bell basis. The fidelity F = F\ is the probability of find- 
ing the pair in the desired Bell state. The shape param- 
eter v = \ (F-^+F^+Fi) measures, e.g., the relative ratio of 
bit error to phase error in the generated pair. We may 
classify various entangled states according to (F,v), and 
use this classification to facilitate the bookkeeping of dif- 
ferent states. For intermediate distances, we only keep 
track of the minimum time (and the associated density 
matrix) for each class of states labeled by (F,v), rather 
than for each possible state. This significantly allevi- 
ates the computational storage requirements of dynamic 
programming. Meanwhile, since we only use a subset of 
states to optimize longer distance pairs, the obtained pro- 
tocol might be a little slower (no more than 10%) than 
the optimal one. 

In our average time approximation, we only keep track 
of the average time for entanglement generation, connec- 
tion and pumping, instead of the full distribution func- 
tion. This underestimates the time for entanglement con- 
nection, because the average time to generate both sub- 
pairs is longer than the maximum of the average individ- 
ual times for two sub-pairs. Fortunately, our compari- 
son between the average time approach and the Monte 
Carlo simulations show that the average times from the 
average-time approximation and from Monte Carlo sim- 
ulations only differ by a factor of about 2 for both the 
BDCZ and the CTSL schemes. In Fig. 6, the time dis- 
tributions from Monte Carlo simulation are plotted. 



